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Abstract 

We investigate the confinement meclianism in tiiree-flavor QCD with imaginary isospin chemical po- 
tentials {fj-u, fJ-d, fJ-s) = {iOT, —i9T,0), using the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
model, where T is temperature. As for three degenerate flavors, the system has Z3 symmetry at 6 = 2it/3 
and hence the Polyakov loop <P vanishes there for small T. As for 2+1 flavors, the symmetry is not pre- 
served for any 9, but (p becomes zero at ^ = ^conf < 2n/3 for small T. The confinement phase defined by 
^ = is realized, even if the system does not have Z3 symmetry exactly. In the 6-T plane, there is a critical 
endpoint of deconfinement transition. The deconfinement crossover at zero chemical potential is a remnant 
of the first-order deconfinement transition at 9 = ^conf- The relation between the non-diagonal element 
Xus of quark number susceptibilities and the deconfinement transition is studied. The present results can be 
checked by lattice QCD simulations directly, since the simulations are free from the sign problem for any 9. 
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I. INTRODUCTION 



Lattice QCD (LQCD) simulations indicate that QCD is in the confinement and chiral symmetry 
breaking phase at low temperature (T) and in the deconfinement and chiral symmetry restoration 
phase at high T. Understanding of the confinement mechanism is, nevertheless, not adequate, 
although the chiral restoration is relatively well understood. The reason mainly comes from the 
fact that there is no exact symmetry for the deconfinement transition and hence the order parameter 
is unknown. In the limit of infinite current quark mass, the Polyakov-loop UJ] is an exact order 
parameter for the deconfinement transition, since Z^v^ symmetry is exact there, where is the 
number of colors. The chiral condensate is, meanwhile, an exact order parameter for the chiral 
restoration in the limit of zero current quark mass. In the real world where u and d quarks have 
small current masses mi = = m^, the chiral condensate is considered to be a good order 
parameter for the chiral restoration, but there is no guarantee that the Polyakov-loop ^ is a good 
order parameter for the deconfinement transition. 

In order to answer this problem, we constructed a gauge theory invariant under the Z^r^ trans- 
formation, that is, a gauge theory with Nf degenerate flavor fermions having flavor-dependent 
fermion boundary conditions flll |3|]. This Z^r^ -symmetric gauge theory is constructed as follow. 
Let us start with iVc-color QCD. The partition function Z in Euclidean spacetime is 

Z = j DqDqDAexp[-So] (1) 



with the action 

' ' (2) 



where qj is the quark field with flavor / and current quark mass mj, D,y = dy—iA^ is the covariant 
derivative with the gauge field A^,, g is the gauge coupling and F^,^ = d^Ai, — d^A^ — i[Afj_, Ay] = 
F^JV^ for the S\]{Nc) generators T'*. The temporal boundary conditions for quarks are 

g/(x,/3 = l/T) = -g;(x,0). (3) 

The boundary conditions are changed into 

qf{x, 13) = - exp {i2T^k/N,)qf{x, 0) (4) 

by the Zat^ transformation [lli|4KZ|] 

q Uq, 

Ay UAyU-^ -i{dyU)U-\ (5) 



where U{x, r) are elements of SU{Nc) with the property U{x, (3) = exp (—i27ik/Nc)U{x, 0) for 
integer k, while the action 5*0 keeps the original form (|2]) since Zat^ symmetry is the center sym- 
metry of the gauge symmetry [i6J. Z^v^ symmetry thus breaks down through the fermion boundary 
condition in QCD. 

Now we consider the SU(A^) gauge theory with degenerate flavor quarks, i.e. N = Nj = N^, 
and assume the following flavor dependent twist boundary conditions (TBC): 

qf{x,(3) = -exp(-i%)g/(x,0) 

= _ exp [-z{9, + 2n{f - 1)/N)]qj{x, 0) (6) 

for flavors / labeled by integers from 1 to A^; here 9i is an arbitrary real number in a range of 
< < 2n. The action Sq with the TBC is invariant under the Zjy^ transformation. In fact, the 
Zjy^ transformation changes / into f — k, but f — k can be relabeled by / since Sq is invariant 
under the relabeling. This is the gauge theory proposed in our previous works Jll Isl and is referred 
to as the Z^v^ -symmetric gauge theory in this paper. . 
When the fermion field g/ is replaced by 

Qf exp {-iefTT)qf (7) 

with Euclidean time r, the action Sq is changed into IM 

SiOf) = j d^x^ QfhuD^ - l^fli + mf)qf + 



.2 n 



(8) 

with the imaginary quark number chemical potential yU/ = iTOf, while the TBC is transformed 
back to the standard one dS). The action Sq with the TBC is thus equivalent to the action S{6f) 
with the standard one ©. 

In the limit of T = 0, the Z^v^ -symmetric gauge theory is identical with QCD with the standard 
boundary condition (|3]). Noting that <P = at T = 0, one can predict that in the Zjv^- symmetric 
gauge theory Z^r^ symmetry is preserved up to some temperature Tc and spontaneously broken 
above T^.. In fact, this behavior is confirmed in Refs. |2, 3|] by imposing the TBC on the Polyakov- 
loop extended Nambu-Jona-Lasinio (PNJL) model js MBTI] that has the same global symmetries 
as QCD at imaginary chemical potentials. In the Za?^- symmetric gauge theory, the quarkyonic 
phase with ^ = and_finite quark-number density appear at small T and large real quark-number 
chemical potential [3J- ^tv^ symmetry is thus essential for emergence of the quarkyonic phase. 



In real QCD, however, s-quark has a current quark mass heavier than mi, where m/ is a 
current quark mass of light quarks. This means that real QCD does not become Zs-symmetric even 
if the TBC is imposed. However, the TBC or its small extension may be important to understand 
the confinement mechanism, since real QCD with 2+1 flavors is not far from QCD with three 
degenerate flavors. As an extension of the TBC, one can consider the ^^-variant TBC 

{eu,ed,es) = {e,-e,o), (9) 

where 9 varies from to tt; this boundary condition is illustrated in Fig. [T] The system with the 
^^-variant TBC is equivalent to the system with imaginary isospin chemical potential 

{fiu,fid,f^s) = {i0T,-ieT,O), (10) 

where /ii = /Uu — /id = 2i9T means the imaginary isospin chemical potential. The ^^-variant TBC 
agrees with the standard boundary condition (|3} when = and the TBC ^ when 9 = 2n/3. 
For the case of small T, as shown later in Sec. Hill (p keeps a real number for any 9 and varies from 
a positive value to a negative one as 9 increases from 0, since charge-conjugation (C) symmetry 
is not spontaneously broken there. This means that (p becomes zero at some value 6'conf of 9. 
Hence the static-quark free energy, — T ln[^], diverges there. Thus the confinement appears, even 
if the system does not have Z^v^ symmetry exactly. This fact means that one can consider the 
confinement by using (p and regard it as an order parameter of the confinement/deconfinement 
transition, even if Zjv^ symmetry is not preserved exactly. Thus 2-1-1 flavor QCD with imaginary 
isospin-chemical potentials (flOl) is an important system to be analyzed. 

As for zero and real fi, meanwhile, an important question is what is a good indicator of QCD 
phase transition in heavy-ion collisions. About 10 years ago, it was proposed that correlations 
and fluctuations (or susceptibilities) of conserved charges may be good signals [l38l-l43ll. If there 



exists a critical endpoint (CEP) of chiral transition [144 



may be a good signal, since it diverges at the CEP [138| 



'ntnf 



S\, the baryon number susceptibility 



4211 . More recently it was suggested 



that third moments of conserved charges can serve as probes of the CEP ll46ll . Furthermore, it 
was proposed in Ref. [|43n with the two-phase description of quark-hadron phase transition that 
the non-diagonal element Xus of quark number susceptibilities may be a good indicator of the 
confinement/deconfinement transition, since it vanishes in the deconfinement phase where inter- 
actions are weak but becomes finite in the hadron phase where different species of quarks are 



confined. The non-diagonal element was recently analyzed by LQCD simulations 1147 . 
PNJL model |3, 3]. 



48|] and the 



In this paper, we analyze QCD at imaginary isospin-chemical potentials (ITOI) for both cases of 
three degenerate flavors and 2+1 flavors, using the PNJL model. As for three degenerate flavors, 
QCD with imaginary isospin-chemical potentials (fTOl ) becomes Zs-symmetric at 6* = 27r/3 and 
hence <1> becomes zero there. As for 2+1 flavor quarks, <P vanishes at 6^ = 6'conf < 2n/3 and ^^conf is 
a function of T. In the 9-T plane, there is a line of <P = for both three degenerate flavors and 2+1 
flavors. On the line, the thermodynamic potential has the same property between the two cases. 
For real QCD with 2+1 flavors, it is well known that the deconfinement transition is crossover at 
zero 9 115 ill . The deconfinement crossover is a remnant of the first-order deconfinement transition 
at = 6'conf- This 9 dependence indicates that there exists a CEP of deconfinement transition in 
the 9-T plane. 

We also investigate the interplay between the deconfinement transition and the non-diagonal el- 
ement Xus of quark number susceptibilities at zero and finite 9. The pseudocritical temperature of 
the deconfinement transition is usually defined by the peak position of the Polyakov-loop suscep- 
tibility. The absolute value \dxus/dT\ has a peak at the pseudocritical temperature, whereas Xus 
does not. This behavior is more conspicuous at the CEP than at zero 9. In the 9-T plane, hence, a 
transition line defined by the peak position of \dxus/dT\ almost coincides with a transition line by 
the peak position of the Polyakov-loop susceptibility. 

This paper is organized as follows. In §2, the three-flavor PNJL model is recapitulated and 9 
dependence of <P is analyzed with the PNJL model. In §3, numerical results are shown. Section 4 
is devoted to a summary. 



Fig. 1: The ^-variant TBC: Location of e f (/ = u, d, s) on a unit circle in the complex plane. 



5 



II. PNJL MODEL 



The PNJL model [|8l- lllL Il3l- l29l l3ll-l37ll is designed to describe the confinement mechanism as 

reproduc( 

imm, 



The model is also successful for the imaginary isospin chemical potential region [|24 
recapitulate the model. 

The three-flavor PNJL Lagrangian is defined in Euclidian spacetime as 

8 



57 



321 



well as the chiral symmetry breaking. In the imaginary fi region, the model can reproduce LQCD 
data [52-62], since the model describes the Roberge-Weiss (RW) periodicity 



36D 



M- We 



a=0 



detgi(l + 75)'?i + h-c. 



+ l(i<P[A],<P*[A],T), 



(11) 



where Dy = dy — iS^^A^, is the Gell-Mann matrices, rh = diag(m„, m^, ?7is) denotes the 
mass matrix and p, = diag{fiu, fJ-d, fJ-s) stands for the chemical potential matrix. Gs and Gb are 
coupling constants of the scalar-type four-quark and the Kobayashi-Maskawa-'t Hooft (KMT) 



interaction 163 



6411 . respectively. The KM' 



Polyakov-loop <P and its conjugate (p* u, 

1 



interaction breaks f/A(l) symmetry explicitly. The 
650 are defined by 
1 



tr,(L), <P* = -tre(L), 



(12) 



with L = exp(iA4/T) and its hermitian conjugate L. In the PNJL model, they are calculated with 
the Polyakov gauge. We take the Polyakov potential of Ref. (l^: 



U = T^ 



a{T) 



<P*<P 



aiT) 



+ h{T) ln(l - 6^^* + 4(^3 ^ ^*3^) _ 3(<j^*)2) 

To\3 



ao + ai 1^— 



T 



(13) 
(14) 



Parameters of U are fitted to LQCD data at finite T in the pure gauge limit. The parameters except 
To are summarized in Table U The Polyakov potential yields the first-order deconfinement phase 
transition at T = Tq in the pure gauge theory [|66[ 1670 . The original value of To is 270 MeV 
determined from the pure gauge LQCD data, but the PNJL model with this value yields a larger 
value of the pseudocritical temperature Tc at zero chemical potential than T^ ~ 160 MeV predicted 



by fuU LQCD ||68l-|7lll. We then rescale To to 195 MeV so as to reproduce T^ ~ 160 MeV MM- 
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02 


^3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential sector determined in Ref. lll3h . All 
parameters are dimensionless. 



The thermodynamic potential (per volume) is obtained by the mean-field approximation as [|26ll 



n 



■2 E 

f=u,d,s 



P 



(27r)3 

+ UM{of)+U{<P,T), 



(15) 



where 



(16) 
(17) 



with Of = {(jfqf), = Ef ± fif and Ef = yp'^ + Mj . The term J^j {J^j) comes from a quark 
(antiquark) loop. For imaginary chemical potential, <P* is the complex conjugate to ^, since Q is 
real. For imaginary isospin chemical potential, Q is invariant under the C transformation <P -ir^ <P*, 
because 



(18) 



where the second transformation is the relabeling of u and d. The three-dimensional cutoff is taken 
for the momentum integration in the vacuum term [|26|. The dynamical quark masses Mj and the 
mesonic potential t/M are defined by 



(19) 
(20) 



Mf = nif — AGsCTf + 2GY)(7f/afii , 

f=u,d,s 

whcrcfy^f',f^f"md f'^f". 

The NJL sector of the PNJL model has six parameters, (Gs, Gd, rriu, md, rris. A). A typical 
set of the parameters is obtained in Ref. [17211 : for example, m; = m„ = = 5.5 MeV and 



rris = 140.7 MeV. The parameter set is fitted to empirical values of t]'- and vr-meson masses and 
TT-meson decay constant at vacuum. We refer to this parameter set as "set R" (realistic parameter 
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set). For theoretical interest, furthermore, we vary the strange quark mass only as = 5.5 and 
600 MeV. We refer to the parameter set with = 5.5 MeV as "set S" (symmetric parameter set) 
and that with = 600 MeV as "set H" (heavy strange quark parameter set). These parameter 
sets are summarized in Table [III 



Set 


mi (MeV) 


^^(MeV) 


yl(MeV) 






R 


5.5 


140.7 


602.3 


1.835 


12.36 


S 


5.5 


5.5 


602.3 


1.835 


12.36 


H 


5.5 


600.0 


602.3 


1.835 


12.36 



TABLE II: Summary of the parameter sets in the NJL sector. 



The resultant thermodynamic potential 12 is a function of (Xu, a^t, a.^, $ and Values of the 
mean fields are determined from the location of the global minimum of in the variable space. 
Since a ,, = Ud is preserved except for the RW phase in which C symmetry is spontaneously 



broken [12211 . we use a = ((t„ + 0"^ + crs)/3 and a' = ag — (cTu + cr(i)/2 = ctu — o'd as order 
parameters, and adopt <Pr = {<P + <P*) /2 and = (^ - (p*) /{2i) instead of <P and The order 
parameter (Pi is C-odd, whereas a, a' and <Pp are C-even. 



The susceptibilities are calculable as [|17|1 



(21) 

with the curvature matrix C defined by 

" dip.dif, ' ^ ^ 

where ipi = a, = cr', (fs = ^4, = ^i- When it is not confusing, we denote the set 
{/Xu, yUrf, fig} by /i and the set {a, a', ^r, <Pi} by (p. The Polyakov-loop susceptibilities, x<Pr<Pr ^^d 
X<Pi<i>i are real for imaginary isospin chemical potential except for the RW-phase. 
The susceptibilities of quark number densities are defined as 

D'Q{T,fi,cp{T,f^)) 

^ ^J^' ' ^ ^ 

where /, /' = u,d,s and the derivation means the partial derivation with respect to fif with 

fixing the other external parameters T and /^/"(t^ ///). We also define the derivative of the suscep- 
tibilities with respect to T as 

Xff',T = • (24) 

8 



Using the Z3 transformation 

18,22,32,3]: 



with an arbitrary integer k, one can see that H has the RW periodicity [|6, 



f^{0^, Od, Os) = n{9u + 2kn/3, 9d + 2kn/3, 9^ + 2k7i/3). 



(26) 



The RW periodicity does not mean that the system is Zs-symmetric, since the external parameters 
6*/ are shifted by the Z3 transformation. An exception is the three degenerate flavor system with 
9 = 2n/3. In fact, the 6*/ are shifted by the Z3 transformation, but the shifted 9f are transformed 
back to the original by the relabeling of flavors, as mentioned in Sec. HI In the confinement phase 
appearing at low T as a consequence of Z3 symmetry, the thermodynamic potential is obtained by 
setting <P = in (fTSi) . The resultant thermodynamic potential includes only three-quark configu- 
rations e~^^f^^. As a result of this property, the flavor symmetry broken by the flavor-dependent 
TBC is recovered in the confinement phase M>\M- 



Finally we consider ^^-dependence of <P at low T. In (fT5l) . i? depends on only through the 
Polyakov-loop potential W and the logarithmic term F = — T ^^(In J-/ + In J^)- For small T, it 
is well satisfied that T <^ Mj and ^ <^ 1. Hence F is approximated into 

F ^ -6T<PN (27) 

with 

N = ^e-^^y, (28) 
/ 

where <P and N are real for any 9, because C symmetry is preserved. For small T, the Polyakov- 
loop potential U has a global minimum at<P = 0. The logarithmic term then moves the minimum 
point to positive (negative) when N is positive (negative). When 9 varies from to 27r/3, N 
changes the sign from plus to minus. Hence the minimum point moves from positive ^ to negative 
<!> as 9 increases from to 27r/3. Eventually the confinement phase with <P = emerges at some 
value 6'conf of 6^. It follows from = that 

^conf ^ arccos f-^e^^^^'^'A , (29) 
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and hence ^conf = 27r/3 for three degenerate flavors (M; = Ms) and 7i/2 for two flavors (Ms = 
oo). This indicates that 6'conf is between 2n/3 and n/2 for 2+1 flavors, whereas 6'conf = 2tt/3 for 
set S. Thus the confinement emerges, even if Z3 symmetry is not preserved. The thermodynamic 
potential (fT5l) in the phase has only three-quark configurations e'^^f^"^ for both of three degenerate 
flavors and 2+1 flavors. The only difference between the two cases is the value of ^^conf • 

ni. NUMERICAL RESULTS 
A. The case of = 

We first analyze the thermodynamics at = through the PNJL calculation with set R. In panel 

(a) of Fig. [2l cr and <P are plotted as a function of T. As T increases, a decreases smoothly, while ^ 
increases continuously. Both the chiral and the deconfinement transition are thus crossover. Panel 

(b) corresponds to T dependence of susceptibilities Xaa, X<Pr<Pr ^^id X'Pi'Pi - Here the pseudocritical 
temperature Tq (Td) of the chiral (deconfinement) transition is defined by the temperature at which 
Xacr (x<Pb,<Pb) becomes maximum; in the present calculation, Tq = 202 MeV and Td = 161 MeV. 
Since ^ = <Pi is zero. This means that <Pi itself is not a good order parameter of the deconfine- 
ment transition. The peak position of x<Pi<Pi does not coincide with that of x<pR<fR' but it should be 
noted that x<Pi'Pi is rapidly changed at T = Td where x<pR<Pn ^ peak. 




T (GsV) T (GsV) 

Fig. 2: T dependence of (a) the chiral condensate a and the Polyakov loop ^ and (b) their susceptibilities 
Xctct . X^R<Pn ^^'^ X^i<Pi at = 0. Set R is taken in the PNJL calculation. The chiral condensate a is 
normalized by the value (Jq at T = 0. X^>b^b ^'^^ X^'j'Pi multiplied by 10 and 100, respectively. 
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B. Thecaseof 6* = 27r/3 



Next we analyze the thermodynamics at ^ = 27r/3 through the PNJL model with set R and set 
S. First we consider set R. In Fig. HJ^a), a and |^| are plotted as a function of T. The Polyakov loop 
<P is real at small T where C symmetry is preserved, but becomes complex at high T where the C 
symmetry is spontaneously broken. The high T region is called the RW phase; further discussion 
will be made in subsection IIII CI for the RW phase. As T increases, <P has a discontinuity at 
T = Td = 188 MeV, whereas a decreases slowly. This indicates that the deconfinement transition 
is the first order. Thus the imaginary isospin chemical potential makes the deconfinement transition 
stronger. The same property is seen in the Nj = 2 case ll24ll . 





Chiral condensate — — 


(a) 


Polyakov- loop 









Chiral condensate — — 




Polyakov-loop 







Fig. 3: T dependence of the chiral condensate a and the absolute value |^| at 6 = 27r/3. The PNJL 
calculation is done with set R in panel (a) and set S in panel (b). The chiral condensate a is normalized by 
the value do at T = 0. 

In the RW phase, the isospin symmetry between u and d is also broken due to the spontaneous 
breaking of C symmetry [22]. Figure Ufa) shows T dependence of af at 9 = 2n/3 for the case of 
set R. As a consequence of the isospin symmetry breaking, one of light quarks becomes heavier, 
while the other becomes lighter. 

Figure ^h) is the same as Fig. IHa), but the PNJL calculation is done with set S. The order 
parameters a and <P have similar T dependence to those in Fig. [3ta). However, <P is zero below 
Td = 194MeV, since Z3 symmetry is exactly preserved. Z3 symmetry is spontaneously broken 
above T^. In this case, the confinement/decconfinement phase transition is govemed by Z3 sym- 
metry. Figure Hfb) is the same as Fig. Ha), but the PNJL calculation is done with set S. In panel 
(b), two of three quarks are always degenerate and become heavier at T > Td, while one of them 
becomes lighter there. The flavor symmetry broken by the TBC is found to be recovered in the 
confinement phase below Td. 
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0.02 - 



0.015 ^ 

0.01 - 

0.005 - 

— 
0.1 



Chiral condensate (light 1) 
Chiral condensate (light 2) 
Chiral condensate (strange) 



(b) Chiral condensate (light) 

Chiral condensate (heavy) 




Fig. 4: T dependence of the chiral condensates cij at ^ = 27r/3. The PNJL calculation is done with set R in 
panel (a) and set S in panel (b). In panel (a), the solid line represents s-quark, while the dashed and dotted 
lines correspond to light quarks. In panel (b), the soUd and dashed lines agree with each other. 

C. Phase diagram in the 6-T plane 

Comparing Fig. |2] with Fig. |3l one can predict that there exists a critical endpoint (CEP) 
of deconfinement transition in the 0-T plane. In fact, the CEP appears at = 6'cep = 0.62 x 
27r/3 = 0.4l7r, since the transition changes the order from crossover to first order there. Figure |5] 
shows T dependence of x#r<?Sr' X<?i#i ^^d Xao- ax9 = 6'cep- The susceptibility X<?'h<z^r diverges at 
T = Tqep = 175 MeV. This indicates that the deconfinement phase transition is second order just 
on the CEP. Since <P^ and a are C-even, they can be correlated with each other. In fact, Xo-a has a 
divergent peak on the CEP, although it has another peak at higher T. Meanwhile, <Pi is C-odd, so 
that it is not divergent on the CEP, although it rapidly decreases there. 





Polyakov-loop real 

'olyakov-loop imaginary — — . 
Chiral condensate — — 


- 


- 



T (GeV) 



Fig. 5: T dependence of Xua, X#r*h ^rid x^\^\ at 6* = 6'cep- The PNJL calculation is done with set R. Xaa 
is multiplied by -1. x#r$r is divided by 10, while x-^i^i is multiplied by 100. 

As mentioned in Sec. [Ill <P becomes zero at some value 6'conf of Q when T is small, and ^^conf is 
27r/3 for three degenerate flavors and less than 27r/3 for 2+1 flavors. Figure |6] shows Q dependence 
of ^ at low T. The PNJL calculation is done for three cases of set R, set S and set H. The value 
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of ^conf is 27r/3 for set S (m, = 5.5 MeV), ~ 0.85 x 27r/3 = 0.567r for set R (m, = 140.7 MeV) 
and ~ 0.75 X 27r/3 = 7r/2 for set H (m^ = 600 MeV). The set-H case agrees with the two-flavor 
case [24]. Strictly speaking, 6'conf depends on T for 2+1 flavors, but does not for three degenerate 
flavors, as shown in (|28T ). 





SetR 




Sets 




SetH - - 


\ 























Fig. 6: 9 dependence of ^ at T = 150MeV. The PNJL calculations are done with three parameter sets of R, 
S and H. 

Figure |7] shows Tdependence of the phase of at 6* = 27r/3. The phase </> is C-odd and hence 
an order parameter of C symmetry. In panel (a) of set R, is tt at small T, and it jumps to about 
±27r/3 at high T. This indicates that ^ is negative at small T and C symmetry is broken at high 
T. Thus the system is in the RW phase at high T. Similar result is seen in panel (b) of set S. Here 
<P is zero at small T and hence is not defined there. At high T, is ±27r/3. 



(b) 

2 - 



Fig. 7: T-dependence of the phase of # at = 27r/3. The PNJL calculation is done with (a) set R and (b) 
set S. Note that = 47r/3 = -27r/3 mod 2-k. 

T dependence of at = 27r/3 shown above is illustrated in Fig. [8] as a movement of <P with 
respect to increasing T in the complex ^ plane. In panel (a) of set R, the RW phase transition as 
the spontaneous breaking of C symmetry is illustrated by two arrows from some negative value. 
The transition is mirror- symmetric about the real ^ axis in virtue of C symmetry. In panel (b) of 
set S, the deconfinement transition as the spontaneous breaking of Z3 symmetry is illustrated by 
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three arrows. The transition is Z3- symmetric and a jump of ^ = to ^ = (p'e^^'^/^ for positive 
<!>' smaller than 1 . 




Fig. 8: Typical movement of at = 27r/3 in (a) set R and (b) set S. The Polyakov loop <P is not allowed 
to move outside the dashed hnes, since the Polyakov-loop potential diverges on the dashed lines. 

Figure |9] shows 9 dependence of (p at high T. The 9 dependence is symmetric with respect to 
the line = tt, so we consider a range of < 9 < n only. In panel (a) of set R, (p is zero at 
< 9 < 0.91 X 27r/3 = O.GItt, but diverges into two solutions (f) ^ ±2n/3 at 9 > 0.61n. The 
system is in the RW phase when 9 > O.Glvr; here one solution is a C image of the other solution. 
Similar result is seen in panel (b) of set S; is zero at < 9 < 27r/3, while ~ ±27r/3 at 
9 > 2n/3. 

The RW phase transition was discovered first in the 9g-T plane 1^], where 9q is the dimension- 
less imaginary quark-number chemical potential defined by the quark-number chemical potential 
fig as /ig = iT9q. The quark-number density is discontinuous at 9g = vr/S mod 2n/3, when T is 
higher than some temperature Trw. This first-order phase transition is now called the RW phase 



transition. It was found in Ref. [1221] that C symmetry is also spontaneously broken in the RW 
phase transition. One can then define the RW phase transition by the spontaneous breaking of C 
symmetry. In the 9g-9-T space, the RW phase appears as a plane of 6'^ = 0, 6*0 < < vr and 
T > Tpfw. Here 9q is a critical value of 9 and depends on T for set R but not for set S; for example, 
9o = 0.93 X 27r/3 at T = 250 MeV for set R and 27r/3 at T > Trw for set S. 

Figure \T0\ shows the phase diagram in the 9-T plane. Three panels correspond to parameter 
sets of R, S and H, respectively. In panel (a) of set R, as mentioned above, there is a CEP of 
deconfinement transition at ^^cep = 0.4l7r and Tcep = 175 MeV. Meanwhile, the chiral transition 
is crossover. Phase diagrams for set S and set H have almost the same structure as that for set 
R, as shown in panels (b) and (c). The dot-dashed line (the left boundary of RW phase) and the 
dot-dot-dashed line (the line of$ = 0) are shifted to the left, when becomes heavy. The shift of 
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Fig. 9: 6'-dependence of the phase (pof^ atT = 250MeV. The PNJL calculation is done with (a) set R and 
(b) set S. 

the dot-dot-dashed line indicates that the confinement at 6* = becomes stronger as rus increases. 
In all the cases, the chiral restoration is weakened at large 6, so the chiral-transition lines are not 
shown there. 







(a) 




/ 






RW 


^ 1 






*<0 


* >0 



(b) 






RW 














* >0 


*<0 



/ 

(c) / 

y 


\ 

RW 


*>0 


'I><0 



Fig. 10: Phase diagram in the Q-T plane. The PNJL calculation is done with (a) set R, (b) set S and (c) set H. 
The solid, dotted, dashed lines represent first-order deconfinement, crossover deconfinement and crossover 
chiral transitions, respectively. The dot-dashed line is the left boundary of the RW phase, while the lower 
boundary of the RW phase is the solid line. The dot-dot-dashed line is the line on which ^ = 0. Below the 
RW region, ^ becomes negative. 

In Fig. [m \^\ is plotted as a function of Q and T. For large \^\ has an abrupt jump as T 
increases. This jump indicates that the deconfinement transition is first order. As Q decreases, 
the order of the deconfinement transition is changed into crossover. The crossover deconfinement 
transition at = is thus a remnant of the first-order deconfinement transition at large Q. 
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Fig. 1 1 : The absolute value of as a function of 9 and T. The PNJL calculation is done with set R. 



D. Susceptibilities at zero and finite 9 

Figure [l2ta) shows T dependence of the non-diagonal element Xus of quark number density 
susceptibilities and its derivative Xus,t with respect to T for the case of 6* = 0. As an interesting 
property, Xus suddenly changes near T-q = 161 MeV, and consequently |Xus,t| has a peak there. 
Comparing Xus in Fig. [I2ta) with x<i>i<Pi in Fig. [2tb) shows that that the two quantities are strongly 
correlated with each other. 

The interplay between Xus and x<Pi<Pj can be understood as follows. The stationary condition 

d 



■f?(T,/i,<^) = (30) 



and its derivatives 

D 



lead to a relation between Xff and x^PiiPj as 



d 

n{T,fi,i^{T,fi)) 



dip 



(31) 



Xff 



DfifDfifi 

^ d'f2{T,ix,^) ^^ d'niT,fx,^) ^ ^^^^ 

dfifd(fi "^'"^^ d^f'difj 

For later convenience, we define dx = d/dX. As for Xus, the second derivative d^,^d^^n vanishes 
at fif = 0, since is /xj-even for each / and hence df^^d^^il is /i/-odd. Similarly, /ij-odd 
quantities d^^-d^Q , 8^^.8^1(2 and d^^.d,i,^Q vanish at /i/ = 0, whereas a /i/-even quantitiy d^^d,p-^n 
does not. These properties lead to 
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Thus Xus is correlated with not x<pR<fR but x<Pi<Pi- The PNJL result shown in Fig. [I2ta) well sim- 
ulates T dependence of Xus calculated with LQCD [1471 14811. although the former underestimates 
the latter for the magnitude of Xus- This underestimation may stem from the fact that the present 
model does not treat baryon degrees of freedom explicitly. 
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Fig. 12: T dependence of quark number susceptibility Xus (solid line) and its derivative (dashed line) with 
respect to T for (a) 6* = and (b) 6 = Oqep- The PNJL calculation is done with set R, and the resultant 
Xus is multiplied by 20. The corresponding LQCD result shown by dots with errorbars in panel (a), 
are not normalized. 

Similar behavior is seen on the CEP, as shown in Fig. H^ b). The relation (|33l) persists for 
finite 9, since C symmetry is preserved there. Therefore |x«s| rapidly decreases at T = Tcep, and 
consequently Xus,t has a peak there. The peak is even more conspicuous for 9 = 9cep than for 
^ = 0. 



IV. SUMMARY 



We have investigated the confinement mechanism in three-flavor QCD with finite imaginary 
isospin chemical potentials (yU^, /id, /is) = {i9T, —i9T,0), using the PNJL model. Two cases 
of three degenerate flavors and 2-1-1 flavors were taken. In the former case, the system has Z3 
symmetry at = 2n/3. The symmetry is not spontaneously broken for low T and hence the 
Polyakov loop ^ is zero there. In the latter case, Z3 symmetry is explicitly broken for any 9, but 
<P becomes zero at 6* = ^^conf smaller than 2tt/3. Thus the confinement phase defined by ^ = 
is realized, even if the system does not have Z3 symmetry exactly. In the confinement phase, the 
static-quark free energy, — T ln[^], diverges. This means that one can consider the confinement by 
using <P and regard it as an order parameter of the confinement/deconfinement transition, even if 
Z3 symmetry is not preserved exactly. 
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The phase diagram is determined in the 9-T plane. There is a CEP of deconfinement transition 
at a finite value 6'cep of 6. This is a good contrast with the CEP of chiral transition at real quark- 
number chemical potential. As another interesting point, there is a line of = for both cases 
of three degenerate flavors and 2+1 flavors. On the line, the confinement phase is realized and 
the thermodynamic potential has only three-quark configurations in which red, green and blue 
quarks are statistically in the same state. This property is independent of whether the system has 
Z3 symmetry or not. 

In the 0-T plane, the confinement transition is crossover on the axis of 6* = 0, while it is first 
order on the line of = 0. This means that the deconfinement crossover at 6* = is a remnant 
of the first-order deconfinement transition on the line of ^ = 0. Hence the distance of the line 
of ^ = from the axis of 6' = shows how strong the confinement property is at 6* = 0. The 
line of ^ = is moved to smaller 6, as nis increases. This means that the confinement property 
is stronger in the two-flavor case than in the 2+1 flavor case. This statement may be consistent 



with LQCD results [133| 



68L 



7l|l . The non-diagonal element Xus of quark number susceptibilities 



is correlated with not X<?r<?r but x<i>i<^i for zero and finite 6. Hence \xus \ does not have a peak near 
the pseudocritical temperature defined by x<2>h<?'r' but rapidly changes there. As a consequence of 
this property, its derivative Xus,t has a peak there. The peak is more conspicuous for finite 6 than 
for zero /i. 

The present results are derived with the PNJL model, but these can be checked by LQCD 
simulations directly, since the simulations are free from the sign problem for imaginary isospin 
chemical potential. 
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